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We report on several methodologically different NNLO determinations of the strong coupling 
constant from event shapes and related observables in e+e- annihilation. The event shape distri- 
butions are analyzed within a combined framework of NNLO and NNLA resummation. We also 
investigate the role of hadronisation corrections, using both Monte Carlo generator predictions 
and analytic models to parametrise non-perturbative power corrections. 
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1. Introduction 

The reaction of e+e^ annihilation into 3 jets has played historically a very prominent role 
for phenomenology and allows a precise determination of the strong coupling constant as, since 
the deviation from two-jet configurations is proportional to it. The phenomenologically interesting 
observables for the purpose of determining the strong coupling constant can be divided into two 
big categories. On one side there are jet rates, which rely on a jet algorithm and are based on the 
number of jets observed in an event, on the other side also the topology of the single events can 
be studied in a systematic fashion by means of so-called event-shape observables. The popularity 
of the latter is mainly due to the fact that they are well suited both for experimental measurement 
and for theoretical description, since many of them are infrared and collinear safe. They describe 
topological properties of hadronic final states by parameterizing the energy-momentum flow of an 
event. This class of observables is also interesting for the study of hadronisation effects. Hadro- 
nisation corrections usually result in a distortion of the event-shape distributions and are thus not 
easily disentangled from partonic predictions. However, the study of mean values and higher mo- 
ments of event shapes allows a clear separation of perturbative and non-perturbative predictions. 
Furthermore, the comparison of hadronization corrections computed by general purpose Monte 
Carlo programs with predictions from analytical models permits to shed some light on the effects 
of hadronization corrections in the determination of a^. 

We have studied jet rates, the distributions and the first five moments of six event-shape ob- 
servables: thrust T (respectively, T = 1 — T), heavy jet mass p, wide and total jet broadening 
B^/v and Bj, C-parameter and the two-to-three-jet transition parameter in the Durham algorithm, 
73. For the definitions of these variables and their historical origin we refer to [|I|] and references 
therein. We will denote the variables collectively as y in the following, such that their two jet limit 
is J — )• 0. Recently, with the appearance of the NNLO results, several extractions of a^ have been 
performed using event-shape observables [^ ^ ^ ^ ^. Their results are summarized in Fig. |l|. In 
the following we present three of them ||3|, ^ ^. 

2. Us from Event- Shape Distributions 

The fixed-order QCD description of event-shape distributions is given by a perturbative ex- 
pansion of the form 



1 d(7 , , _ , , dA , , _ 9 , , dB , , _ , , , dC 

- — -T-{y,Q:^) = as{n)—{y) + at{^) — {y,xy,) + a's{^)— 

c^had ^y ^y ^y ^y 



{y,x^) + i?{at), (2.1) 



where (Xs = ^ and x^ = ^, and where A, B and C are the perturbatively calculated coefficients [^] 
at LO, NLO and NNLO. The distribution is normalised to the total hadronic cross section CJj^^j in 
e^e^ annihilation. The dependence of (2.1) on the collision energy is only through a, and x^ and 



the scale dependence of a^ is determined according to the three-loop running of UsilJ.). 

In order to obtain reliable predictions over the full kinematical range, the perturbative fixed- 



order prediction (2.1) has to be matched with resummation, which is taken into account at next- 
to-leading logarithmic (NLL) accuracy in the In /^-matching scheme. For more details about the 
NLLA-i-NNLO matching we refer to Ref . [S] and references therein. 
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We have used the six event-shape observables listed in Section |I| for our fits. The measure- 
ments we use have been carried out by the ALEPH collaboration [§] at eight different centre-of- 
mass (CM) energies between 91.2 and 206 GeV. The perturbative QCD prediction is corrected for 
hadronisation and resonance decays by means of a transition matrix, which is computed with the 
MC generators PYTHIA [|I0||, HERWIG [||] and ARIADNE |[l|], all tuned to global hadronic ob- 



servables at Mz [|13|]. Corrected measurements of event-shape distributions are compared to the 
theoretical calculation at particle level. For a detailed description of the determination and treat- 
ment of experimental systematic uncertainties we refer to Refs. [0, |2l 



The combined results of six event-shape variables and eight LEP1/LEP2 CM energies is 

a,(Mz) =0.1224 ± 0.0009 (stat) ± 0.0009 (exp) ± 0.0012 (had) ± 0.0035 (theo) . 

For the fitted values of the coupling constant as found from event-shape variables calculated at var- 
ious orders we refer to the figures and tables of |^. The central value of the result is slightly lower 
than the central value of 0. 1228 obtained from a fit using purely fixed-order NNLO predictions [^]. 
Furthermore the dominant theoretical uncertainty on as{Mz), as estimated from scale variations, is 
reduced by 20% compared to NLO-i-NLLA and the scatter among the values of as{Mz) extracted 
from the six different event-shape variables is substantially reduced. However, compared to the 
fit based on purely fixed-order NNLO predictions, the perturbative uncertainty is increased in the 
NNLO-i-NLLA fit. The reason is that in the two-jet region the NLLA-i-NLO and NLLA-i-NNLO 
predictions agree by construction and therefore the renormalisation scale uncertainty is dominated 
by the resummation in this region, which results in a larger overall scale uncertainty in the a, fit. 
Apart from the a^ determination using the standard MC generators mentioned above, we used 



HERWIG++ [jlj version 2.3 together with the MCNLO [jlj] and POWHEG [[Ij] schemes for in- 
vestigating hadronization corrections. From this study it appears that there are two "classes" of 
variables. With standard hadronisation corrections from PYTHIA we obtain as{Mz) values some 
5% higher for the first class, consisting of T , C and Bj, which still suffer from sizable missing 
higher order corrections, than for the second class consisting of the p, Bw and Y^, which have a 
better perturbative stability [|l^]. This means that the PYTHIA hadronisation corrections, applied in 
the ttv fit, might be too small for the first class of variables, resulting in a larger as{Mz) value. For 
further details of our analysis we refer to Ref. [^. 

3. Us from Moments of Event Shapes 

The nth moment of an event-shape observable y is defined by 

{f) = -— y'^dy, (3.1) 

o^had ^0 ay 

where ymax is the kinematically allowed upper limit of the observable. For moments of event 
shapes, one expects the hadronisation corrections to be additive, such that the cross section can be 
divided into a perturbative and a non-perturbative contribution, where the non-perturbative contri- 
bution accounts for hadronisation effects. 

In ref. [Q], the dispersive model derived in Refs. |18, [T^, 20] has been used and extended 



to NNLO to estimate hadronisation corrections to event-shape moments by calculating analytical 
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predictions for power corrections. It introduces only a single new parameter Oq, which can be 
interpreted as the average strong coupling in the non-perturbative region: 

- [ 'dQa,s{Q^) = ao{^i), (3.2) 

IJ-i Jo 

where below the IR cutoff Hj the strong coupling is replaced by an effective coupling. This disper- 
sive model for the strong coupling leads to a shift in the distributions 

where the numerical factor ay depends on the event shape, while P is believed to be universal and 



scales with the CM energy like jJ-i/Q. Insertion of eq. (33) into {3A) leads to 



(/) = / dy{y + a,Pr—-^{y) « / dy{y + ayPr — -^{y) . (3.4) 

J-ayP Otoi ay Jo Giot ay 

From this expression one can extract the non-perturbative predictions for the moments of y. 

The expressions derived in |0] match the dispersive model with the perturbative prediction 
at NNLO QCD. Comparing these expressions with experimental data on event-shape moments, a 
combined determination of the perturbative strong coupling constant a^ and the non-perturbative 
parameter Oq has been performed [Q], based on data from the JADE and OPAL experiments [pT|]. 
The data consist of 18 points at CM energies between 14.0 and 206.6 GeV for the first five moments 



of r, C, F3, p, Bw and Bt, and have been taken from []22[]. For each moment the NLO as well as 
the NNLO prediction was fitted with av(Mz) and Oq as fit parameters, except for the moments of 
73, which have no leading ^ power correction and thus are independent of Oq. 

Compared to previous results at NLO, inclusion of NNLO effects results in a considerably 
improved consistency in the parameters determined from different shape variables, and in a sub- 
stantial reduction of the error on a^. Furthermore the theoretical error on the extraction of as{Mz) 
from p, Yj and Bw is considerably smaller than from T, C and Bj. As mentioned above, the mo- 
ments of the former three shape variables receive moderate NNLO corrections for all n, while the 
NNLO corrections for the latter three are large already for « = 1 and increase with n. Consequently, 
the theoretical description of the moments of p , F3 and B^ displays a higher perturbative stabihty, 
which is reflected in the smaller theoretical uncertainty on a^(Mz) derived from those variables. 

In a second step, we combine the as{Mz) and Oq measurements obtained from different event- 
shape variables. Taking the weighted mean over all values except Bw and Bj, we obtain at NNLO: 

a,(Mz)= 0.1153 ±0.0017(exp)±0.0023(th), ao = 0.5132 ±0.0115(exp)±0.0381(th). 

The moments of Bw and Bt have been excluded here since their theoretical description requires an 
additional contribution to the non-perturbative coefficient P [Q] which is unknown to NNLO. 

The average of a^iMz) is dominated by the measurements based on p and Y3, which have the 
smallest theoretical uncertainties. From NLO to NNLO [Q], the error on as{Mz) is reduced by 
a factor of two. The error on as{Mz) is clearly dominated by the x^ variation, while the largest 



contribution to the error on Oq comes from the uncertainty on the Milan factor ^ [19], which has 
not been improved in the current study. 
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Figure 1: Recent determinations of a^ based NNLO predictions taken from Refs. 



To quantify the difference of the dispersive model to hadronisation coiTections from the legacy 
generators, we analysed the moments of (1 — T) with hadronisation corrections from PYTHIA. As a 
result, we obtained fit results for tts (Mz) which are typically 4% higher than by using the dispersive 
model, with a slightly worse quality of the fit. Comparing perturbative and non-perturbative contri- 
butions at ^/s = Mz, we observed that PYTHIA hadronisation corrections amount to less than half 
the power corrections obtained in the dispersive model, thereby explaining the tendency towards a 
larger value of av(Mz), since the missing numerical magnitude of the power corrections must be 
compensated by a larger perturbative contribution. 



4. a^ from Jet Rates 

Very recently a new determination of a^ was performed using jet rates ^. Theoretical NNLO 



predictions for jet rates [ |23| ] as a function of the jet resolution parameter jcut are compared to 
ALEPH data [^ using the Durham jet algorithm, for which the distance measure is given by 



y,-j,D = 2 min{Ef, EJ) {\ - cos d,-j) /E[ 



2 

vis 



(4.1) 



where E^^^ denotes the energy sum of all particles in the final state. In Ref . [ 23 1 it was shown that 



NNLO predictions for jet rates have only very small hadronization corrections and the theoretical 
error for 10^' < 3'cut < 10^^ drops below the per-cent relative uncertainty. This motivates a ded- 
icated extraction of a^. The corrected ALEPH measurements for the three-jet rate are compared 
to the theoretical calculation at particle level. Values for a^(Mz) are obtained by a least-squares 
fit, performed separately for each jcut value in the range . A nice stability of the result is found up 
to values of ln3'cut ~ —4.5. As final result the value for jcut = 0.02 is taken, which represents an 
optimal compromise between minimal systematic uncertainty and stabihty. The following values 
for as is found: 

a,(Mz) =0.1175 ±0.0020(exp)±0.0015(th), . (4.2) 

Results from LEP2 energies give similar central values but larger statistical uncertainties. Com- 
bining the errors in quadrature yields af:{Mz) = 0.1175 ±0.0025, which is in excellent agreement 
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with the latest world average [|^]. This verifies the expectations that the three-jet rate is an ex- 
cellent observable for this kind of analysis, thanks to the good behaviour of its perturbative and 
non-perturbative contributions over a sizable range of jet-resolution parameters. 

5. Acknowledgments 

This research was supported in part by the Swiss National Science Foundation (SNF) under 
contracts PP0022- 118864 and 200020-126691. 

References 



[1 

[2 
[3 

[4; 

[5 
[6 

[7 

[8 

[9 

[10 

[11 

[12 

[13 
[14 

[15: 
[16 
[17 
[18 

[19 
[20 
[21 
[22 
[23 
[24; 



R. W. L. Jones et al, JHEP 0312 (2003) 007 [arXiv:hep-ph/0312016]. 

G. Dissertori et al, JHEP 0802 (2008) 040 [arXiv:0712.0327]; 
S. Bethke et al, Eur. Phys. J. C 64 (2009) 351 [arXiv:0810.1389]. 

G. Dissertori et al, JHEP 0908 (2009) 036 [arXiv:0906.3436]. 

T. Gehrmann, M. Jaquier and G. Luisoni, Eur. Phys. J. C 67 (2010) 57 [arXiv:091 1.2422]. 

G. Dissertori etal, Phys. Rev. Lett. 104 (2010) 072002 [arXiv:0910.4283]. 

T. Becher and M. D. Schwartz, JHEP 0807 (2008) 034 [arXiv:0803.0342]; 

Y. T Chien and M. D. Schwartz, [ai'Xiv: 1005. 1644]; 

R. Abbate, M. Fickinger, A. H. Hoang, V. Mateu and I. W. Stewart, [arXiv: 1006.3080]. 

A. Gehrmann-De Ridder et al, JHEP 0712 (2007) 094 [arXiv:071 1.471 1]. 

T. Gehrmann, G. Luisoni and H. Stenzel, Phys. Lett. B 664 (2008) 265 [arXiv:0803.0695]. 

A. Heister et al [ALEPH Collaboration], Eur. Phys. J. C 35 (2004) 457. 

T. Sjostrand et al, Comput. Phys. Commun. 135 (2001) 238 [hep-ph/0010017]. 

G. Corcella et al JHEP 0101 (2001) 010 [hep-ph/001 1363]. 

L. Lonnblad, Comput. Phys. Commun. 71 (1992) 15. 

R. Barate et al [ALEPH Collaboration], Phys. Rep. 294 (1998) 1. 

O. Latunde-Dada, JHEP 0711 (2007) 040 [arXiv:0708.4390]; 

O. Latunde-Dada, S. Gieseke and B. Webber, JHEP 0702 (2007) 051 [hep-ph/06 12281]. 

S. Frixione and B. R. Webber, JHEP 0206 (2002) 029 [hep-ph/0204244] . 

P Nason, JHEP 0411 (2004) 040 [hep-ph/0409146]. 

A. Gehrmann-De Ridder et al, JHEP 0905 (2009) 106 [arXiv:0903.4658]. 

Y. L. Dokshitzer et al, Nucl. Phys. B 469 (1996) 93 [hep-ph/9512336]; 

Y L. Dokshitzer and B. R. Webber, Phys. Lett. B 404 (1997) 321 [hep-ph/9704298]. 

Y L. Dokshitzer et al, JHEP 9805 (1998) 003 [hep-ph/9802381]. 

Y L. Dokshitzer et al, Eur. Phys. J. direct C 1 (1999) 3 [hep-ph/98 12487]. 

C. Pahl et al, Eur. Phys. J. C 60 (2009) 181, (E) C 62 (2009) 451 [arXiv:0810.2933]. 
C. J. Pahl, Ph.D. thesis, MPI Munich 2007. 

A. Gehrmann-De Ridder et al, Phys. Rev. Lett. 100 (2008) 172001 [arXiv:0802.0813]. 
S. Bethke, Eur. Phys. J. C 64 (2009) 689 [arXiv:0908.1135]. 



